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Abstract 

The stability of GPS time and frequency transfer is limited by the fact that 
GPS signals travel through the ionosphere. In high precision geodetic time transfer 
(i.e. based on precise modeling of code and carrier phase GPS data), the so-called 
ionosphere-free combination of the code and carrier phase measurements made on 
the two frequencies is used to remove the first-order ionospheric effect. In this paper, 
we investigate the impact of residual second- and third-order ionospheric effects on 
geodetic time transfer solutions i.e. remote atomic clock comparisons based on GPS 
measurements, using the ATOMIUM software developed at the Royal Observatory 
of Belgium (ROB). The impact of third-order ionospheric effects was shown to be 
negligible, while for second-order effects, the tests performed on different time links 
and at different epochs show a small impact of the order of some picoseconds, on a 
quiet day, and up to more than 10 picoseconds in case of high ionospheric activity. 
The geomagnetic storm of the 30th October 2003 is used to illustrate how space 
weather products are relevant to understand perturbations in geodetic time and 
frequency transfer. 
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1 Introduction 



Ti me and freq uency transfer (TFT) using GPS (Global Positioning Sys- 
tem, ( Leick . 2OO4I )) or equivalently GNSS (Global Navigation Satellite System) 



satellites consists in comparing two remote atomic clocks to the reference time 
scale of the GPS system (or to another post-processed time scale based on a 
GPS or GNSS network). From the differences between these comparisons, one 
gets the synchronization error between the two remote clocks and its time 
evolution. 

GPS TFT is widely used within the time community, for example for the re- 
alization of TAI (Temps Atomique International), the basis of the legal time 
UTC (Universal Time Coordinated), computed by the Bureau International 
des Poids et Mesures (BIPM). TFT is characterized by its very good reso- 
lution (1 observation point/30s or possibly 1 point per second) and a high 
precision and frequency stability thanks to the carrier phases (uncertainty uA 
of about 0.1 ns). The present uncertainty in GPS equipment calibration is 5 
ns (uncertainty uB -systematic, hence calibration errors- in the BIPM circular 
T). 

The solar activity varies according to the famous 11 -year solar cycle, also 
called sunspot cycle. An indicator of the solar activity is the number and in- 
tensity of solar flare events, as large flares are less frequent than smaller ones 
and as the frequency of occurrence varies. Solar flares are violent explosions in 
the Sun's atmosphere th at can release as much e nergy as 6 • 10^^ Joules (more 



details can be found in (IBerghmans et al.l . l2005l )). The X-rays and UV radi- 



ation emitted by the strongest solar flares can affect the Earth's ionosphere, 
modifying the density and the distribution of the electrons within it. 
Some solar flares give rise to Coronal Mass Ejections (CMEs), i.e. ejections of 
plasma (mainly electrons and protons) from the solar corona, carrying a mag- 
netic field and travelling at speeds from about 20 km/s to 2700 km/s, with an 
average speed of about 500 km/s. Some CMEs reach the Earth as an Inter- 
planetary CME (ICME) perturbing the Earth's magnetosphere and hence the 
Earth's ionosphere. Consequently, solar flares and associated CMEs strongly 
influence our terrestrial environment which has an impact on GPS/GNSS sig- 
nals. 

Since ionospheric influence on electromagnetic waves is frequency depen- 
dent and since GPS signals are broadcasted in two different frequencies, iono- 
spheric effects are commonly removed through a given combination (named 
ionosphere- free) of the signals in the two frequencies /i and f2- However, it 
is well known that this combination removes only first-order perturbations, 
which correspond to about 99.9% of the total perturbation. The present study 
aims at evaluating the impac t of the remain i ng pa rt, c oncentrating on second- 



and third-order effects. While iFritsche et al.l (120051 ) and iHernandez-Paj ares et al 
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( 120071 ) investigated the higher-order ionospheric impact on GPS receiver po- 
sition, GPS satelhte clock or GPS satelhte position estimates, in the present 
paper, we focus on the impact on precise time and frequency transfer using 
GPS signals. 

Second- and third-or der ionospheric t e rms a re therefore implemented in the 



software ATOMIUM flDefraigne et all 120081 ). developed at the Royal Obser- 



vatory of Belgium. ATOMIUM is based on a least-square analysis of dual- 
frequency carrier phase and code measurements and is able to provide clock 
solutions in Precise-Point-Positioning (PPP) as well as in single-difference 
(also called Common- View, CV) mode. 



The present paper is organized as follows. The next section recalls the 
principles of GPS TFT and the ionosphere-free analysis in Precise-Point- 
Positioning or Common- View mode. In Section 3, the ionosphere-free analysis, 
as implemented in the ATOMIUM software, is reviewed. In Section 4, the se- 
lected method used to implement higher-order ionospheric corrections in the 
ATOMIUM ionosphere-free analysis is described. Our corresponding results 
are presented in Section 5, in terms of ionospheric delays of second and third 
orders compared to first-order ionospheric effects, and then in terms of the 
impact of higher-order ionospheric delays in the receiver clock solution com- 
puted with ATOMIUM. In Section 6, we put our results in perspective with 
some investigations on related solar flare events and K-index considerations. 
Conclusions are finally presented in Section 7. 

Numerical values in the equations of this paper are provided in SI units. 



2 GPS time and frequency transfer 



High precision geodetic GPS receivers can lock their internal oscillator on 
an external frequency, given by a stable atomic clock. The GPS measurements 
are then based on the clock frequency. Using post-processed satellite orbits and 
satellite clock products computed by the International GNSS Service (ICS), 
one can deduce the synchronization error between the external clock and either 
the GPS time scale or the reference time scale of the ICS, named IGST. The 
use of GPS measurements made in two stations p and q gives then access to 
the synchronization error between the two atomic clocks in these stations p 
and q. 

For a station p or similarly q, the GPS measurements, relative to observed 
satellite i, on the signal code Pk and phase Lk, at frequency k (1 for /i = 
1575.42 MHz or 2 for /2 = 1227.6 MHz) with corresponding wavelength A^, 
can be written in length units as 



3 



k=1.2 



pi + cAtp + cAt' + Ar; + (epj; + (+71^ + 2 ■ 72^ + 3 • 13^ 



[Lk=i,2)l = p; + cAtp + cAt' + Ar; + N;\k + (e^j; + (-7U - 12, 



p 

(1) 



following iBassiri fc Hajjl (119931 ) who considered separately the different orders 
of the ionosphere impact on GPS code and phase measurements. There, 
is the geometric distance i-p ; Atp is the station clock synchronization error; 
Ar* is the satellite clock synchronization error; Ar^ is the tropospheric path 
delay for path i-p; Ilk (oc 72fc (oc 1//|) and 73fc (oc 1//^) are the first-, 

second- and third-order ionospheric delays on frequency k; are the phase 

ambiguities; (e^Pj.)^ and (^Lj.)^ are the error terms in code P and phase L, 
respectively, containing noise such as unmodeled multipath, hardware delays 
(biases from the electronic of the satellite or the station receiver) on the prop- 
agation of the modulation/carrier of the signal. In their most general form, 
the different orders of ionospheric delays (71 ,72, 73) are expressed as inte- 
grals along the true path of the signal (including a bending, which is function 
of the signal frequency, in the dispersive ionosphere), in terms of the signal 
frequency, of the local electron densi ty in the ionosphere a nd of the local ge- 
omagnetic field along the trajectory (jBassiri &: Hajjl . 1 19931 ). 
When a dual frequency GPS receiver is available at station p, the so-called 
i onosphere- fr ee combination [k = 3) is used. This combination is defined as 
flLeickl . l2004 ): 



^3 = 
L3 = 



n 



if? - /I) 



Pi- 



Li 



/I 



if? -fir 



fi 



(2) 



with /i and /2 the two GPS carrier frequencies. This combination removes, 
from the GPS signal, the first-order ionospheric effect, 71, since the latter 
is proportional to the inverse of the square frequency. The corresponding 
ionosphere-free observation equations therefore do not contain any first-order 
ionospheric ter m, but contain new factors for second- and third-order iono- 
spheric effects (IFritsche et al.l . l2005l ) with respect to previous Equation 1: 



[P^)l = pl + cAtp + cAr* + Arl + {ep^ + (+2 •723 - 3-73 



(L3); = pl + cAt, + cAt' + Ar; + N;X, + (5^3); + (-723 + IhT, (3) 



Note that about 99.9% ( Hernandez- Pa j ares et al. . 20081 ) of ionospheric per- 
turbations are removed with 71 in the so-called ionosphere-free combination. 
Note also that while the first order has the same magnitude on GPS phase and 
code measurements (but with opposite sign), the impact of second- and third- 
order effects is larger on code than on phase observations (twice for 72, three 
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times for 13). The GPS observation equations given by 3 are directly used in 
Precise Point Positioning. The clock solution obtained is the synchronization 
error between the receiver clock and the GPS or IGS Time scale. 



For Common- View analysis, one uses the single differences between simul- 
taneous observations of a same satellite i in two remote stations p and q in 
order to determine directly the synchronization error between the two remote 
clocks: the observation equations for receivers p and q with satellite i are 
subtracted. This single difference removes the satellite clock bias in the GPS 
signal, assuming that the nominal times of observation of the satellite by the 
two stations are the same. When forming ionosphere-free combinations, the 
single-difference code and carrier phase equations are: 



(A);, = + cAtp, + Ar;^ + {ep,l^ + (+2 • 723 - 3 • 13,);^ 

(Lslq = + cAtpq + Ar;, + N;^X, + (5^3);^ + {-12^ + /Ss);, (4) 

where, for any quantity X, 



Now, all the terms in above Equations 1, 3 or 4 can be estimated via an 
inversion procedure using some a priori precise satellite orbits and satellite 
clock products. This finally provides the solution for either Atp in PPP, i.e. 
the clock synchronization error between the atomic clock connected to the GPS 
receiver and the GPS or IGS Time scale at each epoch, or Atpq in Common 
View, i.e. the synchronization error between the remote clocks connected to 
two GPS receivers. In parallel, the station position and tropospheric zenith 
delays are estimated as a by-product. 



3 The ATOMIUM software 



The present study on ionosph eric higher-order pert urbations in TFT is 



based on the ATOMIUM software flDefraigne et al.l . l2008l ). which uses a weighted 



least-square approach with ionosphere-free combinations of dual-frequency 
GPS code (P3) and carrier phase (L3) observati ons. ATOMIUM was ini tially 



developed to perform GPS PPP, as described in (iKouba &: Herouxl . l200ll ). and 
later adapted to single differences, or Common View (CV), of GPS code and 
carrier phase observations. In the paragraphs below, we describe ATOMIUM, 
following the diagram presented in Figure 1. First, GPS ionosphere-free code 
and phase combinations are constructed according to Equations 2 from Li, Pi, 
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L2, P2 observations read in RINEX observation files. By default, the ATOM- 
lUM software uses as a priori the IGS products (liGS data and products 



20091 ) ■ IGS satellite clocks (tabulated with a 5-minute interval) are used to 
obtain Ar* at the same sampling rate as provided. IGS satellite orbits (tabu- 
lated with a 15-minute interval) are used to estimate Pp (or p^^ ) via a 12-points 
Neville interpolation of the satellite position every 5 minutes. The station po- 
sition is corrected for its time variations due t o degree 2 and 3 solid E arth 
tides as recommended by the lERS conventions (iMcCarthv fc Petitl. l2003l. Ch . 
7) and for ocean loading according to the FES2004 model (ILyard et al.l . ]2004h . 
The relative (if GPS observations were made before GPS week 1400) or ab- 
solute (if after) elevation (no azimuth) and nadir dependant corrections for 
receiver and satellite antenna phase center variations are read from the IGS 
atx file available at ( IGS atA 20091 ). Prior to the least-square inversion, the 
computed geometric distance is removed from both phase and code ionosphere- 
free combinations. Those are also corrected for a relativistic (periodic only) 
delay and an hydrostatic tropospheric delay. Tropospheric delays are modeled 
as the sum of a hydrostatic and a wet delay, resulting from the product of a 
given mapping function and of the corresponding hydrostatic or wet zenith 
path d elay (zpd). For the h ydrostatic part, we use the Saastam oien a prior i 
model (ISaastamoinenl . Il972l ) and the dry Niell mapping function flNiell Il996h . 
For the wet part, we use the wet Niell mapping function flNielill996h while the 
wet zpd is estimated as one point every 2 hour, and modeled by linear interpo- 
lation between the se points. Carrie r phase measurements are further corrected 
for phase windup (IWu et al.l . 119931 ) taking into account satellite attitude and 
eclipse events. The implementation of additional higher-order ionospheric cor- 
rections on phase and code is done at this level, as corrections applied on the 
code and phase measurements (see Section 4). 



The least-square analysis used in ATOMIUM is detailed in flDefraigne et al. 



20081 ). As output, ATOMIUM provides the station p (or relative p-q) position 
for the whole day, the receiver clock p (or relative p-q) synchronization error 
every 5 minutes, tropospheric wet zenith path delays p (and q) at a given rate 
(2 hours in our case). 



4 Corrections for ionospheric delays 



First-, second- and third-order ionospheric terms are function of the Slant 
Total Electron Content (STEC), which is the integrated electron density in- 
side a cylinder column of unit base area between Earth ground and satellite 
altitude, along the satellite i- station p direction. STEC is function not only 
of the satellite elevation or station position, but also of the time of the day, 
of the time of the year, of the solar cycle, and of the particular ionospheric 
conditions. As the II term contains 99.9% of the ionospheric perturbations on 
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the GNSS signal, it can be used to estimate STEC. The second- and third- 
order terms are then computed using this estimated STEC. The STEC in II 
can be determined using the geometry-free combinations of the r neasu rements 
made on the two GPS frequencies, defined as P4 and L4 (ILeickl . |2004| ): 



P4 = -Pl+P2 

L4 = Li — L2 (5) 

The quantities (P4L4) only contain a given combination of ionospheric delays 
on the signals of frequency /i and /2, plus some constant terms associated 
with the differential hardware delays (between /i and /2) in the satellite and 
in the receiver, and the phase ambiguities. 

In the following, for a practical and efficient implementation of higher-order 
ionospheric terms in ATOMIUM, we work in the no-bending approximation, 
meaning that the trajectory considered to compute those ionospheric correc- 
tions on GPS observations is a straight line from satellite i to station p (and 



q). Ind eed, from analytica l estimations ba s ed on results of iHoque fc Jakowski 



(l2008l . formula (31)) and IJakowski et al.l (Il994j . formulas (12-13)), we quan- 



tified the bending effects in the dispersive ionosphere on the ionosphere-free 
combinations of GPS measurements. These effects are due to the excess path 
length for a curved trajectory and to the difference in effective STEC for the 
two GPS frequencies. Under extreme conditions, i.e. a highly ionized iono- 
sphere together with very low satellite elevations, ionospheric bending effects 
might reach the level of the third-order ionosphere correction, which turns out 
to be in the TFT noise as we will show in the following. 



We fu rthermore assume the ionosphere Single Layer Model (SLM) (IBassiri fc Hajj 



19931 ). reducing the whole ionosphere layer through which the GPS signal trav- 
els to a single effective sheet at a given height with the equivalent electron 
content (see Figure 3). 

The above assumptions allow to reduce the general integrals for the ionospheric 
delays to the easily implemented equations given in the following subsections. 



4-1 First- order ionospheric delays 



The first-order ionospheric effect, is given by (IBassiri &: HajjI . Il993l ) 



Ilk = aU- STEC (6) 
with the factor for GPS frequencies 1 and 2 being 
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40.3 

/l,2 



ali,2 = +-?2- (7) 



This implies that the corresponding factors for the ionosphere-free {k = 3) or 
geometry-free [k = 4) combinations are 



al4 



As announced here above, for each pair of code or phase measurements (on 
/i and /2), the geometry- free combination can be used to compute the STEC, 
which is needed in higher-order ionospheric corrections. Neglecting the 12 and 
13 contributions induc ing errors in estimated STEC of the order of 0.1 TECU 




at the most, one gets (iHernandez-Pajares et al.l . 120071 ): 



{STECY 



aid 



— ((L^) — (-P4) )arc without — c • DCBp — c ■ DC 



^ '' ' cycle slips 



(9) 



In the above formula. Pi — P2 Differential Code Biases (DCB) are assumed 
constant during a day, and we read them from the CODE (ICS Analysis 
Center) lONEX files; <> means taking the average. 

Alternatively, STEC can be computed u sing PiP? codes th at have first been 



smoothed with the corresponding phase (IDach et al.l . 120071 ) 
1 



{STECY 



{(P2); - (Pi);} - c ■ DCBp - c ■ DCB' 

^ ^ with phase 



(10) 



This leads to similar results as those obtained from Equation 9 with respect 
to the same DCB product. 

Within ATOMIUM, STEC is computed, using formula 9, for each satellite- 
station pair with a sampling rate of 5 minutes. Figure 2 illustrates STEC 
computed via ATOMIUM for stations BRUS, i.e. Brussels, at latitude 50°28' 
and longitude 4°12', TLSE, i.e. Toulouse at latitude 43°20' and longitude 
1°17', ONSA i.e. Onsala at latitude 57°14' and longitude 11°33'. Figure 2a 
features a quiet-ionosphere day (around a minimum of solar activity), 11th 
March 2007, with its normal local-noon diurnal maximum of STEC. Figure 
2b shows the ionospheric storm of October 30, 2003 and the STEC sensitivity 
to perturbed ionospheric conditions as induced by solar activity (see Section 
6). Since ionospheric effects in GPS are function of STEC, we understand that 
ionosphere- induced errors will increase in the next few years due to the increas- 
ing solar activity associated with the ascending phase of the 24th sunspot cycle 
(maximum forecast around 2011-2012 depending on the models). 
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4-2 Second-order ionospheric delays 



Whereas the magnitude of II for a given frequency depends solely on STEC 
and is always positive, the magnitude and sign of 12 depend on the geomag- 
netic field B values, on STEC and on the i-p signal direction via the angle 
between B and the Line of Sight (LOS), Ob-t.os (Figure 3 ). We used the 



following integrated formula (iHernandez-Pa j ares et all l2007f ) 



I2k = a2k- Bjpp ■ cos 9b-los ■ STEC (11) 
with frequency factors 



7527 • c 



7527 ■ c 

a2i,2 = 73 — (13) 

Jl,2 

STEC is obtained from L4P4 (Equation 9) and Bjpp is con iputed using the ac - 



curate International Geomagnetic Reference (ICR) model (iTsyganenkd . l2005l ) . 
as the latter allows to reduce errors in 12 up to 60% with respect to a dipolar 
model (Hernandez- Pa j ares et al., 2007, 2008). 

4.3 Third-order ionospheric delays 



In the third-order ionospheric contribution, the magnetic field term can be 
safely neglected at s ub-millimeter error level, leading to the simple formula 



fiFritsche et al.l . l2005l ) 



I3k = aSk ■ STEC (14) 

with frequency factors also being functions of the electron distribution in the 
ionosphere: 



2437 • iV^ax • V 

a3i,2 = 74 (15) 

/l,2 



"33 = 3,^2. f2 (16) 
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where the shape factor r] is taken around 0.66 and the peak electron density 
along the signal propagation path, N^ax-, can be determined by a linear in- 
terpolation between a typical ionospheric situation and a solar maximum one 



( iFritsche et al.l . l2005l . formula 14 corrected according to a private communi- 



cation from M. Fritsche, Novembe r 2008), using, as interpo lation coefficients 



the numerical values suggested by lBrunner fc M. Gd (Il99ll ): 



^max = ^\'olV.^,^ ■ {yTEC - 4.55 ■ lO'') + 20 • 10^^ (17) 



[(4.55 - 1.38) ■ 10 



The Vertical TEC (VTEC), which is TEC along a vertical trajectory below 
the satellite, is taken as the projection, via the Modified Single Layer Model 
ionosphere mapping function, of {STEC)p from Equa tion 9 with_aM,s'r/ if 



0.9782, REarth = 6371 • 10^ m, H = 506.7 • 10^ m as in flDach et all . 120071 ): 

STEC = fusLuiz) ■ VTEC (: 



fMSLhliz) = l/i 



REarth ^ ^ 



\ V REarth + H 



■ C0S2 {aMSLM ■ z) (19) 



Finally, using the above equations for the first-, second- and third-order iono- 
spheric effects on GPS signal propagation, one finds the orders of magnitude 
of their impact on the code and carrier phase measurements given in Table 1. 



5 Quantifying ionospheric effects on TFT 



The 12 and /3 corrections computed according the procedure described 
above were applied to the ionosphere-free combinations P3 and L3 used in 
ATOMIUM. The present section shows some preliminary results: estimated 
second- and third-order delays on GPS signals (and on combinations of their 
measurements) and the impact of these on the time and frequency transfer 
solutions. 



5.1 Ionospheric delays 



The first results concern the ionospheric delays as computed with ATOM- 
IUM according to the models detailed in previous section. 
Firstly, recall that the Total Electron Content of the ionosphere reaches its 



10 



maximal value at local noon, on a normal day (Figure 2). Since II, 12 and 13 
are proportional to Slant TEC, the amplitude of ionospheric perturbations in 
GPS signal reflect this daily variation of TEC. 

Secondly, the diurnal TEC maximum is function of the station latitude, so the 
ionospheric delays in GPS signals follow accordingly. 

And finally, for any observed satellite, as the ionospheric thickness crossed by 
the signal is proportional to the inverse of the sine of the satellite elevation, 
the STEC during one satellite track, as well as the ionospheric delays, takes 
the shape of a concave curve. 

Figures 6, 7 and 8 illustrate /I2, 12^ and 13^ respectively on a quiet (left) 
versus an ionosphere-stormy day, the ionospheric storm of October 30, 2003 
(right). The selected station in this illustration is Onsala (ONSA). Each point 
of the curves corresponds to an ionospheric correction on a GPS phase mea- 
surement for an observed satellite. The first-order ionospheric perturbations 
in L2 can reach about 100 nanoseconds during the storm (Figures 6 a and b), 
while it is less than 50 nanoseconds in normal times. II in Li is slightly smaller 
according to factor fi/ff- The II effect is removed from the ionosphere- free 
combination. 

The second-order ionospheric perturbation in the ionosphere-free combination 
(Figure 7 a) is about 3 to 4 orders of magnitude smaller than the first order 
in L2. 123 can reach about 15 picoseconds during the storm, about the double 
of its maximum value during a quiet day (Figure 7b compared to 7a). We also 
recall that, in the ionosphere- free combination, the second-order ionospheric 
perturbation, 12^, affects twice more the codes than the phases, as seen in 
Equation 3. 

Figure 8 illustrates the third-order ionospheric perturbation which is again an 
order of magnitude smaller than the second order. Here the effect of the storm 
is also clear, as the third-order effect in the ionosphere-free combination can 
reach about 2 picoseconds during the storm (Figure 8b), while its maximum 
value on a non-stormy day is about 0.14 picoseconds (Figure 8a). Again, the 
contribution of 13^ is three times more important for codes than for phases, as 
seen in Equation 3, but remains negligible with respect to the present precision 
of GPS time and frequency transfer. 



5.2 Ionospheric impact on receiver clock estimates from a analysis 



Table 1 and the results presented in the above paragraph illustrate the 
need to take second-order ionospheric corrections into account in P and L 
measurements for TFT. However, to be coherent, in addition to the 12 (and 
/3) correction(s) on GPS code and phase data, we should also use satellite 
orbit and clock products computed with 12 (and /3) correction(s) in order 
to estimate the impact of the ionosphere on station clock syn chronization er- 
rors via ATOMIUM. Indeed, iHernandez-Pajares et al.l (120071 ) estimated that 
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second-order ionospheric effects in satellite clocks were the largest and could 
be more than 1 centimeter (i.e. 30 picoseconds); the same authors men- 
tioned that the second-order ionospheric effects on the satellite position are 
of the order of several miUimeters only, and consist in a global southward 
shift of the constellation. Current IGS products do not take 12 nor 73 into 
account, this is why we present here the impact of our ionospheric corrections 
on clock solutions via ATOMIUM in Common- View mode (Figures 9 and 10), 
as the satellite clock is eliminated in CV. We choose the link BRUS-ONSA, 
i.e. Brussels- Onsala (Sweden), the day of an ionospheric storm, October 30, 
2003, and used GPS observations with a sateUite elevation cutoff of 5 degrees. 
Figure 9 shows the effect of applying the 12 corrections on GPS P3L3 analysis. 
We see an effect up to more than 10 picoseconds during the ionospheric storm 
on the hnk BRUS-ONSA. 

The 73 effect shown in Figure 10 is at the present noise level of GPS observa- 
tions. 

Consequently, the residual ionospheric errors in P3L3 (when P3L3 is not cor- 
rected for higher-order effects) are mainly due to the contribution of 123. A 
12 delay of 15 picoseconds peak to peak during the storm (Figure 7b) for a 
given station A induces a variation with the corresponding differential 723^ - 
723s amplitude in CV frequency transfer with station B, as the shape of the 
curve is determined by the GPS phases for which the 72 correction is applied 
with a factor 1. Furthermore, 72 induces twice as much an offset on the ab- 
solute time synchronization error (Figure 9), as the calibration of the curve 
is determined by the code data for which the 72 correction is applied with a 
factor 2 (Equation 3). However, this of course is still well below the present 
calibration capabilities of GPS equipment. 

Note that the results presented here correspond to the time link BRUS-ONSA. 
It is therefore the differential ionospheric effect between those two stations that 
matters for the clock solution in Common- View mode. The impact of 72 on a 
clock solution in PPP could therefore be higher and induce larger effects on 
intercontinental time links. This will be investigated in further studies, when 
consistent satellite orbits and clock products computed with 72 (and 73) will 
be available. 



6 Space weather 

Solar activity was initially monitored thanks to ground observations (sunspot 
group evolution). Since recent decades, the Sun is observed via sateUites. This 

opened the field to space weather studies, characterizing the conditions of the 
Sun, of the space between the Sun and Earth, and on the Earth. 
The Geostationary Operational Environmental Satellite (GOES, orbiting the 
Earth at 35790 km) measures solar flares as X-rays from 100 to 800 picometers 
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and classifies them as A, B, C, M or X according to their peak flux (in watts 
per square meter) on a logarithmic scale. Each class has a peak flux ten times 
greater than the preceding one, with X ones of the order of 10~^ W/m^. There 
is an additional linear scale from 1 to 9 inside each class. 

Furthermore, the magnetic fleld of the Earth is characterized (according 
to Earth latitude) by several indexes. The K index quantifles disturbances 
in the horizontal component of the Earth's magnetic fleld with an integer 
in the range 0-9. K index is derived from the maximum fluctuations of the 
Earth geomagnetic fleld horizontal component observed on a magnetome- 
ter during a three-hour interval. The conversion table from maximum fluc- 
tuation (in nanoTesla) to -ftT-index, varies frorn observatory to observatory 



(INOAA Space Weather Prediction Centerl . 120091 ). Hence, the official planetary 



Kp index is derived by calculating a weighted average of i^-indices from a net- 
work of geomagnetic observatories. Values of K from 5 and higher indicate a 
geomagnetic storm which may hamper GPS/GNSS TFT as we shall see. 

The previous solar maximum was around 2001. The year 2003 was thus 
an active period for the Sun, compared to year 2007. For the sake of com- 
parison, we selected a quiet day, March 11 2007, and the very active day of 
30th October 2003. Around our quiet test-day, no X and no M flares were 
observed; only four B-type X-ray events were detected during March 2007. 
On the other side, the solar activity and Earth geomagnetic conditions were 
at exceptionally high levels around the 30th October 2003, as can be seen 
from Figure 4, which lists the solar flare events according to their intensity, 
together with the Earth geomagnetic index K from Wingst Observatory {Kw 
in the text) in Germany (at geographic latitude 53.74° and longitude 9.07°, 
close to Onsala's latitude). Those events were mostly due to two large solar 
active regions named NOAA0486 and NOAA0484, which produced numerous 
solar M-flares and several X-flares. Further insight for Figure 4 can be gained 
when rea ding the corr esponding SIDC (Solar Influences Data Center) weekly 



bulletins (jSIDQ l2009l ) on solar activity for weeks 148 and 149 of 2003, corre- 
sponding to the range of our plots. We focus on the X-flares accompanied by 
CME. 

On October 28, a XI 7. 2 (thus extremely strong) solar flare occurred from the 
solar active region NOAA0486 peaking at 11:10 UT. It was accompanied by a 
CME directed towards the Earth with an estimated plane-of-the-sky speed of 
about 2125 km/s, flrst detected at 10:54 UT in the LASCO C2 fleld of view 
of the SOHO spacecraft (located at about 1.5 • 10® km from the Earth). As 
the CME was extremely fast, the shock arrived to the Earth around 06:00 UT 
on October 29. It produced a severe magnetic storm with K^r index reaching 
9. The intensity of the storm decreased slightly to Kw = 7 during October 
29 as the arriving magnetic cloud was of the north-south type. The main por- 
tion of the negative the north-south interplanetary magnetic fleld component 
Bz arrived in the trailing part of the cloud at the end of October 29. Bz was 
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strongly negative during around 8 hours producing the second peak in the Kw 
index which reached 9 again during 18:00 - 24:00 UT. The index settled 
down to minor storm conditions {K^ = 5) on October 30. 
On October 29, the solar active region NOAA0486 produced an XlO.l flare 
peaking at 20:49 UT, as observed by GOES. It was associated with a CME 
observed by SOHO/LASCO C2 at 20:54 UT that developed to a full halo 
CME directed towards the Earth with an estimated plane-of-the-sky speed of 
about 1950 km/s. The electromagnetic shock of the CME was registered by 
ACE/MAG around 16:00 UT on October 30. This time the magnetic cloud 
was of south-north type, so the severe geomagnetic storm started right after 
the arrival of the shock: the Kw index reached 9 again and stayed at that 
level during two 3-hour intervals (18:00 - 24:00). The geomagnetic storm fi- 
nally ended on October 31 - November 1 when the Kw index dropped to 4. 
From October 30 to November 5 included, several X- and M-class flares oc- 
curred in solar active regions NOAA0484, NOAA0486 and NOAA0488. 
On November 2, an X8.3 flare was observed by GOES, peaking at 17:25 UT, 
from the solar active region NOAA0486. A halo CME with estimated plane- 
of-the-sky speed of about 2100 - 2200 km/s followed. The shock of the arriving 
CME was recorded in the solar wind on November 4, at 05:53UT. The associ- 
ated interplanetary magnetic fleld pointed southward between 07:00 and 9:30 
UT. This triggered a geomagnetic storm {Kw=G-7), but only for a limited 
duration: thereafter, the geomagnetic field remained quiet to unsettled. 
On November 4, came from NOAA0486 an extreme X17 fiare (peaking at 
19:53UT), later estimated to have reached an X28 peak flux that saturated 
the GOES detectors, together with a full halo CME. The shock, correspond- 
ing to this CME which arrived sideways on November 6 at 19:37 UT, was 
relatively weak and only led to a minor storm that lasted until November 7, 
00:00 UT. The magnctosphcrc then remained quiet to unsettled. 
By contrast, during the end of the week, the Sun had a very low activity. 
However, an existing large solar coronal hole gave rise to minor geomagnetic 
storm {Kiv — 6) episodes on November 9. 

We now come back to our STEC at station Onsala on the 30th October 
2003 (Figure 5), estimated using the software ATOMIUM and GPS data. We 
understand that the huge solar flare event XI 7. 2 of October 28 2003, that was 
associated with a CME, through the Earth geomagnetic storm it triggered on 
the 29th October 2003, is responsible for the abnormally high level of STEC at 
Onsala (and Brussels as well) at early hours of October 30 2003. This explains 
the correspondingly high first-, second- and third-order ionospheric perturba- 
tions computed in Figures 6, 7 and 8, in the first hours of October 30 2003. A 
small corresponding impact can be seen on the ATOMIUM clock solution in 
Figure 9. 

The major Earth geomagnetic storm that occurred on the 30th October 2003, 
however, is most probably due to the solar fiare event XlO.l of October 29 
2003, associated with a CME. This time, the strong signature present in the 
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STEC estimated with ATOMIUM for Onsala (Figure 5) (and Brussels as well) 
and computed corresponding first-, second- and third-order ionospheric pertur- 
bations on GPS signals (Figures 6, 7 and 8) is clearly visible in the estimated 
Common- View clock solution of ATOMIUM (Figure 9). This emphasizes the 
need for higher-order ionospheric corrections on GPS signals in precise time 
and frequency transfer, as we need to prepare for the next solar maximum. 
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7 Conclusions 



This study aimed at quantifying the impact of second- and third-order 
ionospheric delays on precise time and frequency transfer using GPS signals. 

We used the ATOMIUM software, in which we implemented higher-order 
ionospheric contributions (second and third orders) in the so-called ionosphere- 
free combination of GPS codes and phases. We then compared these iono- 
spheric contributions with the first-order ionospheric effect on the GPS dual 
frequency signal, which is cancelled in ionosphere-free combination. It was 
shown that the first-order ionospheric delay of several tens of nanoseconds 
on an ionosphere-quiet day, is doubled in case of ionospheric storms. Though 
second-order delays in the ionosphere-free combination are about 3 to 4 or- 
ders of magnitude smaller than the first-order delays, they can reach about 15 
picoseconds on a stormy day, which is significant when performing geodetic 
time and frequency transfer with very stable clocks. Third-order delays in the 
ionosphere-free combination are yet an order of magnitude smaller, and are 
at the level of present noise of GPS observations. The impact of those higher- 
order delays on ionosphere-free time and frequency transfer clock sohitions was 
estimated for the time link BRUS-ONSA. It reaches more than 10 picoseconds 
during the ionospheric storm of October 30 2003. Finally, we illustrated the 
correlations between space weather data (strong solar events, disturbances of 
the Earth geomagnetic field) and time/frequency transfer performed with GPS 
signals. 
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Fig. 1. ATOMIUM software diagram. 
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Fig. 2. STEC computed when observing GPS satellite prn 18 with GPS receivers 
located at Brussels (BRUS), Onsala (ONSA) and Toulouse (TLSE) on an iono- 
sphere-quiet day (11th March 2007), Figure a (left), or on a stormy day (30th 
October 2003), Figure b (right). STEC was computed with the ATOMIUM soft- 
ware, using Equation (9) and DCBs read from CODE lONEX files. 1 TECU= lO^'' 
e~/m^. 
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Fig. 3. The second-order ionospheric effect is not only a function of STEC, but it 
is also function of the angle between the Line Of Sight (LOS) and the geomagnetic 
field B at the Ionosphere Piercing Point (IPP), and of the magnitude of B at IPP. 
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Fig. 4. Time series of the noticeable solar events (X and M X-ray flares) and geo- 
magnetic K index from Wingst Observatory {Kw in the text). The notation used 
is the following. Labels i and j are linked to the classification of flares. M flares are 
represented by a line peaking between 10 and 20 (e.g. M2.5 by a line from to 12.5), 
while X flares are represented between 20 and 30, at their respective peak-intensity 
time. Note that X flares bigger than X9.9 saturate the scale and are represented by 
a line peaking at 30. 
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Fig. 5. STEC time series computed with the ATOMIUM software for station Onsala 
(ONSA), according to Equation 9, from 26th October to 9th November 2003. 
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Fig. 6. First-order ionospheric delay in GPS frequency 2, for station Onsala on an 
ionosphere-quiet day, 11th March 2007 (Figure a, left), versus on an ionosphere-s- 
tormy day, 30th October 2003 (Figure b, right). 




Fig. 7. Second-order ionospheric delay in GPS so-called ionosphere- free combination, 
for station Onsala on an ionosphere-quiet day, 11th March 2007 (Figure a, left), 
versus on an ionosphere-stormy day, 30th October 2003 (Figure b, right). 




Fig. 8. Third-order ionospheric delay in GPS so-called ionosphere-free combination, 
for station Onsala on an ionosphere-quiet day, 11th March 2007 (Figure a, left), 
versus on an ionosphere-stormy day, 30th October 2003 (Figure b, right). 
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Up to -10 ps peak to peak 
during iono storm 



Difference bet. estimated station clock synchro error with 12 (without 13) 

versus without 12, 13 (seconds). ^% 

The IGS products are used in both cases 
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Fig. 9. Effect of taking second-order ionospheric effects, or not, into account in tlie 
L3P3 GPS measurements for the Brussels-Onsala hnk, on the ionosphere-stormy 
day 30th October 2003. The difference is taken between two ATOMIUM estimated 
station clock solutions, both using IGS products. 
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Fig. 10. Effect of taking third-order ionospheric effect, or not, into account in the 
L3P3 GPS measurements for the Brussels-Onsala link, on the ionosphere-stormy 
day 30th October 2003. The difference is taken between two ATOMIUM estimated 
station clock solutions, both using IGS products. 
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Table 1 

Orders of magnitude of ionospheric effects II, 12 and 13 on GPS phase measure- 
ments (for codes, see converting factor in code measurement Equations 1, 3, Section 

1) 



Ionospheric effect 


Delay in L1L2 


Delay in L3 


Relevance 


(absolute value) 


per 100 TECU 


per 100 TECU 


fHernandez-Paiares et al.. 2008) 


11 


~ 30ns - 100 ns 





99.9% of 1123 


12 


~ - 130 ps 


~ - 45 ps 


90% of 123 


13 


~ — 3 ps 


~ - 2 ps 
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